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Abstract 

A positive semi-definite (PSD) tensor which is not a sum-of-squares (SOS) tensor is 
called a PSD non-SOS (PNS) tensor. Is there a fourth order four dimensional PNS 
Hankel tensor? Until now, this question is still an open problem. Its answer has both 
theoretical and practical meanings. We assume that the generating vector v of the 
Hankel tensor A is symmetric. Under this assumption, we may fix the fifth element 
U 4 of v at 1. We show that there are two surfaces Mq and Nq with the elements 
Vi-, vq, v±,V 3 , V 5 of v as variables, such that Mq > Nq, A is SOS if and only if vq > Mq, 
and A is PSD if and only if vq > No, where vq is the first element of v. If Mq = Nq for 
a point P = (v 2 , vq, v\, vq, u. 5 ) t , then there are no fourth order four dimensional PNS 
Hankel tensors with symmetric generating vectors for such V2, vq, v\, vq, V5. Then, we 
call such a point P PNS-free. We show that a 45-degree planar closed convex cone, a 
segment, a ray and an additional point are PNS-free. Numerical tests check various 
grid points, and find that they are also PNS-free. 

Key words: Hankel tensor, generating vector, sum of squares, positive semi-definiteness, 
PNS-free. 
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1 Introduction 


In 1888, young Hilbert [Sj proved that for homogeneous polynomials, only in the following 
three cases, a positive semi-definite (PSD) polynomial definitely is a sum-of-squares (SOS) 
polynomial: 1) m — 2; 2) n — 2; 3) m — 4 and n = 3, where m is the degree of the polynomial 
and n is the number of variables. Hilbert proved that in all the other possible combinations 
of n and even m, there are PSD non-SOS (PNS) homogeneous polynomials. The most well- 
known PNS homogeneous polynomial is the Motzkin polynomial [TS] with rn — 6 and n = 3. 
Other examples of PNS homogeneous polynomials was found in [1, {3, 4j, 22]. 

A homogeneous polynomial is uniquely corresponding to a symmetric tensor [20]. For 
a symmetric tensor, m is its order and n is its dimension. One class of symmetric tensors 
are Hankel tensors. Hankel tensors arise from signal processing and some other applications 
uummm- ^ pu, two classes of PSD Hankel tensors were identified. They are 
even order strong Hankel tensors and even order complete Hankel tensors. It was proved in 
pH] that complete Hankel tensors are strong Hankel tensors, and even order strong Hankel 
tensors are SOS tensors. It was also shown there that there are SOS Hankel tensors and 
PSD Hankel tensors, which are not strong Hankel tensors. Thus, a question was raised in 
pH] : Are all PSD Hankel tensors SOS tensors? If there are no PSD non-SOS Hankel tensors, 
the problem for determining a given even order Hankel tensor is PSD or not can be solved 
by solving a semi-definite linear programming problem pa mum 

We may call the problem raised by the above question as the Hilbert-Hankel problem, 
as in a certain sense, it is the Hilbert problem with the Hankel constraint. 

According to Hilbert [8], [22], one case with low values of m and n, in which there are 
PNS homogeneous polynomials, is that rn — 6 and n = 3. In Pi. the Hilbert-Hankel 
problem with order six and dimension three was studied. Four special cases were analyzed. 
Thousands of random examples were checked. No PNS Hankel tensors of order six and 
dimension three were found in jX3]. Theoretically, it is still an open problem whether there 
are PNS Hankel tensors of order six and dimension three or not. 

According to Hilbert [8] [22], another case with low values of m and n, in which there 
are PNS homogeneous polynomials, is that m — n — 4. In this paper, we consider Hankel 
tensors of order four and dimension four. 

Let v = (t>o, Ui,..., U12) 1 " G 9ft 13 . A fourth order four dimensional Hankel tensor 
A = (a ili2 i 3 i 4 ) is defined by 

^122^34 ^Zi+22+^3+^4 - 4? 

for 21,22,23,24 = 1,2,3,4. The corresponding vector v that defines the Hankel tensor A is 
called the generating vector of A. For x = (aq, £2, £3, aq) T G 9ft 4 , a Hankel tensor A 
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uniquely defines a Hankel polynomial 


/(x) = A x® 4 = 


4 4 

E \ 

ft ... . rp . rp , rp . rp . - \ n i . . . rp . rp . rp . rp . 

' x ' l 2 dj 1 ' 3 ’ x " l '£ / J ^21+22+23+24—4^21^22 *^24 • 

H^25^3^4 = l *1)^2^3^4=1 


( 1 ) 


If /(x) > 0 for all x G 5+ the Hankel tensor *4, is called positive semi-definite (PSD). 
If /(x) can be represented as a sum of squares of quadratic homogeneous polynomials, the 
Hankel tensor A is called sum-of-squares (SOS). Clearly, A is PSD if it is SOS. 

I 11 the next section, we present some necessary conditions for the positive semi-definiteness 
of fourth order four dimensional Hankel tensors. 

We may see that the role of Vj is symmetric in /(x). In Section 3, we assume that 


Vj = Vn-j 


( 2 ) 


for j = 0,..., 5. Under this assumption, by the results of Section 2, if A is PSD, we have 
v 0 = V 12 > 0 and u 4 = v$ > 0. Moreover, if u 4 = u 8 = 0 and A is PSD, A is SOS. Thus, we 
may only consider the case that v 4 = v 8 > 0. Since A is PSD or SOS or PNS if and only if 
aA is PSD or SOS or PNS respectively, where a is an arbitrary positive number, we may 
simply assume that 

v 4 = v 8 = 1. (3) 

Next, we show that there is a function rj(v 5 ,u 6 ) such that i](v 5 ,v 6 ) < 1 if A is PSD. We 
propose that there are two functions M 0 (v 2 , ^ 6 , ui, U 3 , U 5 ) > N 0 (v 2 ,V 6 ,vi,v 3 ,vs), defined for 
r l{ v bi v &) < U suc h that *4 is SOS if and only if v 0 > M 0 , and A is PSD if and only 
if r > 0 > +)• If Mq = Nq for some U 2 , Ue, V \, U 3 , W 5 , then there are no fourth order four 
dimensional PNS Hankel tensors for such V 2 ,v&, v±, V 3 ,v§ under the symmetric assumption 
([2]). We call such a point P = (u 2 , v 6 , Vi, v 3 , u 5 ) T G 9ft 5 a PNS-free point of fourth order 
four dimensional Hankel tensors, or simply a PNS-free point. We call the set of points in 
+, satisfying ))(v 5 . iq) < 1, the effective domain of fourth order four dimensional Hankel 
tensors, or simply the effective domain, and denote it by S. We show that if all the points 
in S are PNS-free, then there are no fourth order four dimensional PNS Hankel tensors with 
symmetric generating vectors. 

In Section 4, we show that a point P in S is PNS-free if there is a value M, such 
that when vq = M, / 0 (x) = /(x) has an SOS decomposition, and /o(x) = 0 for x = 
(x 4 , X 2 , % 3 , x 4 ) t G 3 ? 4 with x\ + x\ 7 ^ 0. We call such a value M, such an SOS decomposi¬ 
tion of /o(x), and such a vector x the critical value, the critical SOS decomposition 
and the critical minimizer of A at P. respectively. Then, we show that the segment 
L = {(v 2 ,v 6 ,v ll v 3 ,v 5 ) T = (l,l,f,t,f) T : t G [—1,1]} is PNS-free. We conjecture that this 
segment is the minimizer set of both M 0 and N 0 . Then, we show that the 45-degree pla¬ 
nar closed convex cone C = {(u 2 , vq, v±, v 3l vz,) T = (a,b, 0,0,0) T : a > b > 1}, the ray 
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Figure 1: The segment L , the planar closed convex cone C, the ray R and the point A. 

R = {( v 2 , v&, Vi, v 3 , n 5 ) T = (a,0,0,0,0) T : a < 0} and the point A = (1, 0, 0, 0, 0) T are also 
PNS-free. We illustrate L, (7, R and A in Figure |T| 

In Section 5, numerical tests check various grid points, and find that Mq = Nq there. 
Thus, they are also PNS-free. Therefore, numerical tests indicate that there are no fourth 
order four dimensional PNS Hankel tensors with symmetric generating vectors. 

Some final remarks are made in Section 6. 


2 Fourth Order Four Dimensional Hankel Tensors 

We write out ([I]) explicitly in terms of the coordinates of its generating vector v: 

/(x) = Vqx\ + Av\x\x2 + V2(Ax\xz + 6X4X3) + V^{Ax\x\ + Ax\x^ + 12 x\x 2 X%) 

+ V4 (x'2 + 6x^X3 + 12x1X3X3 + 12X3X2X4) + n 5 (4X3X3 + 12x4X2X3 + 12x1X3X4 + 12X3X3X4) 

+ u 6 ( 4 xiX 3 + 4X3X4 + 6x3X4 + 6x3X3 + 24x1X2X3X4) ( 4 ) 

+ U7(4X2X3 + 12X3X3X4 + I2X1X3X4 + I2X1X2X4) + Vg (X3 + 6X3X4 + I2X2X3X4 + I2X1X3X4) 

+ V 9 (4X3X4 + 4X1X4 + I2X2X3X4) + Uio(4x 2 X4 + 6X3X4) + 4^11X3X4 + U12X4. 


The following theorem gives some necessary conditions for fourth order four dimensional 
Hankel tensors being PSD. Particularly, we note that four key elements of its generating 
vector Uo,u 4 ,u 8 ,Ui 2 must be nonnegative. 
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Theorem 1 Suppose that A = (a^; 3 i 4 ) is a Hankel tensor generated by its generating 
vector v = (uo,ui, • • • A 12 ) T G 9ft 13 . If A is a PSD (or positive definite, or SOS, or strong) 
Hankel tensor, then we have 

Vi > 0, (5) 

for 2 = 0 ,4, 8 , 12 , 

Vi + 6 vj +2 + ^+4 > 4|i>j + i + Vi + s\, (6) 

fori = 0,4,8, 

Vi + 6 u i+4 + v i+8 > 41 u i+2 + v i+6 1, (7) 

for 2 = 0,4, and 

v 0 + 6v 6 + vi2>4:\v 3 + v 9 \. (8) 

Proof Let e^ be the 7th column of a 4-by-4 identity matrix, for k = 1, 2, 3,4. Substituting 

x = e k to (QJ for k = 1 , 2 , 3,4, by /(e fc ) > 0 , we have (J5]) for z = 0 ,4, 8 , 12 . 

Substituting x = e*, + e^+i to (QJ for k = 1 , 2 , 3, by /(e^ + e^+i) > 0 , we have 

Vi + 4v i+ i + 6v i+2 + 4u i+3 + v i+4 > 0, 

for i = 0,4, 8 . Substituting x = e*, — e^+i to (QJ for k = 1 , 2 , 3, by /(e^ — e^+i) > 0 , we have 

Vi - Av i+ i + 6v i+2 - 4u i+3 + v i+4 > 0, 

for 2 = 0,4, 8 . Combining these two inequalities, we have (| 6 ]) for 2 = 0,4, 8 . 

Similarly, by /(e*, + ejfc +2 ) > 0 and /(e*, — efc +2 ) > 0 for k — 1, 2, we have (J7J) for i = 0,4. 

By /(e 1 + e 4 ) > 0 and /(e 1 — e 4 ) > 0 , we have (JHJ) - The theorem is proved. □ 

Whereafter, we say that a PSD Henkel tensor is SOS if there is a key element of its 
generating vector Uo,u 4 , Vs,V \2 vanishes. Before we show this, the following lemma is useful. 

Lemma 1 If a polynomial in one variable is always nonnegative: 

p(t ) = aot 2k+1 + a\P k + • • • + u 2 fc+i > 0, V t £ 37. 


Then a 0 = 0. 

Proof If a 0 > 0, we let t —> —00 and get p(t) — > — 00 , which contradicts that p(t) is 
nonnegative. 

If ao < 0, we let t —> +00 and get pit. ) —* — 00 , which also contradicts that p(t) is 
nonnegative. 

Hence, there must be a 0 = 0. □ 

Theorem 2 Suppose the fourth order four dimensional Hankel tensor A is PSD and its 
generating vector is v. If v qUi 2 = 0, then Vj = 0, for j = 1,..., 11, and A is SOS. 
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Proof Without loss of generality, we assume that Vq = 0. 

To prove V\ = 0, we take xi = (t, 1,0,0) T . Then, the homogeneous polynomial (|4|) 
reduces to 

/(x i) = 4vi t 3 + 6 v 2 t 2 + 4 v 3 t + v 4 . 

From Lemma [U we have v\ = 0 since /(xi) is nonnegative. Similarly, we can prove v 2 = 
v 3 = 0 if we take x 2 = (t, 0 , 1 , 0 ) T and x 3 = (t, 0 , 0 , 1 ) T respectively. 

From Theorem Q] we know v 4 > 0. When we take x 4 = (f 2 , t, — 0) T , the homogeneous 

polynomial (S} reduces to 

/(x 4 ) = —( 2 ^ — 2 )u 4 f 4 + (9(f 3 ). 

Let f —> oo. Since /(x 4 ) is always nonnegative, we have u 4 < 0. Hence, there must be 
v 4 = 0. 

If we take x 5 = (t 3 ,0,f, 1) T , the homogeneous polynomial ([4j) is 

/(x 5 ) = 12 v 5 t 7 + 0(t 6 ). 

From Lemma [T] we have v§ = 0 since /(x 5 ) is nonnegative. 

We take x 6 = (t, 0,1,0) T . Then, the homogeneous polynomial (jJJ) is 


/(X 6 ) = 4 : V 6 t + V 8 . 

From Lemma [U we have v§ = 0 since /(xg) is nonnegative. Similarly, we can prove V 7 = 0 
when we take x 7 = ( 0 , f, 1 , 0 ) T . 

We take x 8 = (t 4 , 0, t, 1) T . Then we have 

f (x 8 ) = 12 v 8 t 5 + 0(t 4 ). 

From Lemma [T] we have v 8 = 0 since the polynomial /(x 8 ) is nonnegative. 

We could prove v 9 = 0, v 40 = 0 and i>n = 0 if we takes x 9 = (t, 0, 0,1) T , x 10 = (0, t, 0,1) T 
and Xu = ( 0 , 0 , t, 1 ) T , respectively. 

Finally, since Vo — v 4 — ■ ■ ■ — — 0 , we have 

/(x) = v 12 x\. 

By Theorem [H we get v \ 2 > 0. Hence, the Hankel tensor A is obviously SOS. □ 

Theorem 3 Suppose the fourth order four dimensional Hankel tensor A is PSD and its 
generating vector is v. If v 4 v 8 = 0, then v 3 = 0 for j — 1,2,, 11, and A is SOS. 
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Proof By symmetry, we only need to prove this theorem under the condition u 4 = 0. 

If we take xi = (l,t, 0 , 0 ) T , the homogeneous polynomial (J4J reduces to 

/(x i) = 4 v 3 t 3 + 6 v 2 t 2 + 4vi t + v 0 . 

From Lemma [lj we have v 3 = 0 since /(xi) is nonnegative. Similarly, we can prove v 3 = 
v 6 = 0 if we take x 2 = ( 0 , t, 1 , 0 ) T and x 3 = ( 0 , t, 0 , 1 ) T respectively. 

To prove v-j = 0, we take x 4 = (0,t 2 ,t, 1) T . Then, the homogeneous polynomial (IIJ) 
reduces to 

f (x 4 ) = 16 v 7 t 5 + 0(t 4 ). 

From Lemma [lj we have v~j = 0 since /(x 4 ) is nonnegative. 

From Theorem [H we know v% > 0. When we take x 5 = (0, —t 2 ,t, 1) T , the homogeneous 
polynomial (J4j) reduces to 

/(x 5 ) = -5 v 8 t 4 + 0(t 3 ). 

Let t —>• oo. Since /(x, 5 ) is always nonnegative, we have v$ < 0. Hence, there must be 
v 8 = 0 . 

If we take x 6 = (0, 0, t, 1) T , the homogeneous polynomial (J3|) is 

/(x 6 ) = 4 v g t 3 + 0(t 2 ). 

From Lemma HI we have vg = 0 since /(xg) is nonnegative. Similarly, we could prove r^o = 0 
and Un = 0 if we takes x 7 = ( 0 , t, 0 , 1 ) T and x 8 = ( 0 , 0 , t, 1 ) T , respectively. 

The prove of V\ — 0 and v 2 = 0 could be similarly obtained if we take x 9 = (1, t , 0, 0) T 
and xio = ( 1 , 0 , t, 0 ) T respectively. 

Finally, since v 3 =0 for j = 1,..., 11, we have 

/(x) = v 0 x 4 + v 12 x% 

By Theorem [U we get vq > 0 and v\ 2 > 0. Hence, the Hankel tensor A is obviously SOS. □ 

3 Symmetric Generating Vectors 


Now, we make assumptions ((2]) and Q. At the beginning, we consider a mini problem 
which is the Hankel polynomial with X\ = x 4 = 0. This problem helps us to analyze the 
effective domain of two important surfaces M 0 and N (} . 
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3.1 Function 77 


We consider a two variable quartic polynomial 

9{yi, 2/2) = ayf + 4 :f 3 y :i { y 2 + ^y\y\ + 4 fiy x y\ + 

Its PSD property is completely characterized by the following theorem. 


Theorem 4 The quartic polynomial 3(31, y 2) is PSD if and only if 

[ 4| /31 — 3y ifj <\ P \, 

a>>?(/3,7):= 37 - y V - 8/3= _ lol 

[- 2 - */7> 1/51- 

Proof First, if 2/2) is PSD, from <7(1, — 1 ) > 0 and 3(1, 1 ) > 0 , we have a > 4 |/ 5 | — 3 y. 
Thus, in any case, 77(73,7) > 4 |/ 3 | — 3 y. 

Second, suppose that a > 4 |/ 3 | — 3 y. If 7 < 0 , we get 

3(?/i, 2/2) = (a — 4|/3| + 37)(2/1 + 2/2) + 4|/5|(2 /i + 2/2) 2 (2/1 - 2/12/2 + 2/2) - 3 7(2/? - 2/f) 2 > 0. 


If 0 < 7 < |/ 3 |, we rewrite 27(2/1,2/2) as follows 


3 ( 2 /i, 32) = (a - 4 |/ 3 | + 37)(2/1 + 2/|) + (31 + 32) 2 [( 4 |/ 3 | - 3 y)( 3 2 + 2 /|) - ( 4 |/ 3 | - 67)2/132] ■ 


Since (4|/3| - 67 / - 4(4|/3| - 3y ) 2 = —48|/3|(|/3| - 7 ) < 0, it yi elds that 3 ( 31 , 32 ) > 0. 
Finally, we consider the case 7 > \/3\. Let a = — > q 7 } ierij we } iave 

/ 2/3 x 2 

3(3i, 32) = (« - a)(yi + 32) + a 3? + — 3i32 + yl 

\ a 


Obviously, if a > d, 3(31,32) is SOS and PSD. 

Next, we show that y\ + ^-yvy-i + y\ = 0 has nonzero real roots. For the convenience, we 
denote t — ^ and prove that t 2 + 2 'S-t + 1 = 0 has real roots. It is easy to see that t = 0 is 
not its root. Since 7 > |/ 3 |, we have 


\P\ m 2 |/ 3 |( 3 7 + V 9 7 2 - 8 / 3 2 ) ^ 8 |/ 3 | 7 ^ , 

a 37 - a/ 97 2 - 8 / 3 2 8 / 3 2 8 / 3 2 


Hence, |/ 3 | > a. The discriminant of the quadratic in t is 

OO 2 _4 = 4+^>0. 

\ a ) cr 


Therefore, there are nonzero (31,32) such that 3(31,32) = (a — a)(y\ + 7/2). Obviously, if 
3(31,32) is PSD, we have ct > d. Thus, we say ?/(/?, 7) = d if 7 > |/ 3 |. □ 

Then we have another necessary condition for a fourth order four dimensional Hankel 
tensor A to be PSD under assumptions (J2J) and d 3 J). 












Corollary 1 Under assumptions (U|) and (E|), if A is PSD, then 7 ](v 5 ,vq) < 1. 

Proof Let X\ = X 4 = 0, X 2 = 2/1 and X 3 = 1 / 2 ■ By Theorem HI we have the conclusion. □ 

3.2 Surfaces Mq and No 

We now introduce the key idea of this paper, to establish two surface M 0 and N 0 , in the 
following theorem. 

Theorem 5 Suppose that assumptions (0|) and (E|) hold. Then, there are two functions 
M 0 (v 2 , v 6 , ni, v 3 , v 5 ) > N 0 (v 2 , v 6 , v 1} v 3 , v 5 ) > 0 defined for 


V(v5,v 6 )<l, (9) 

such that A is SOS if and only if v 0 > M 0 (v 2 , v 6 , v 3 , v 3 , v 5 ), and A is PSD if and only if 
Vo > N 0 (v 2 , vq, Vi, v 3 , V5). If for all v 3 and v$ satisfying (EP, we have M 0 (v 2 ,v o,vi,v 3 ,v 3 ) = 
No(v 2 ,V 6 ,vi,v 3 ,V 5 ), then there are no fourth order four dimensional PNS Hankel tensors 
under assumption (U|). 

Proof Using assumptions ([2D and (ED, we rewrite (J3D as 

/(x) = v 0 (xf + x\) + v^fxz + X 3 ) + /i(x) + / 2 (x), 


where 

fi (x) = rj(v 5 , v & ){x\ + xj) + 4 ^ 5 ( 2 : 23:2 + 2 : 2 X 3 ) + 6 v e xlxl 

and 

= 1 - r](v 5 ,v 6 ). 

Then h 4 > 0 by ©• By Theorem [H /i(x) is PSD. Since /i(x) has only two variables, it is 
also SOS by Hilbert p, -221. 

We now consider terms in / 2 (x). Each monomial in / 2 (x) has at least one factor as a 
power of X\ or X4. We may order the monomials of /^(x). For example, consider 12 ^ 52 : 13 : 2 X 3 . 
Assume that it is ordered as the kth monomial of / 2 (x). Then by the arithmetic-geometric 
inequality, we may see that 

-12v 5 X!X 2 xl < 3|v 5 | ^xf + e k xf 2 + 2e k xi}j , 

where e k is a small positive number. We may let e k be small enough such that the sum of 
the coefficients for 2:2 on the right hand side of the above inequality for all possible k is less 
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than v±. By symmetry, the sum of the coefficients for x\ on the right hand side of the above 
inequality for all possible k is less than V 4 . We see that 

12v 5 X!X 2 xl + 3|n 5 | (x\ + + 2e k xl 

\ e k 

is a PSD diagonal minus tail form. By [ 6 ], it is SOS. Thus, as long as Vq is big enough, when 
(JSJ is satisfied, /(x) is SOS. From this, we see that M 0 and N 0 exist, such that they are 
defined as long as AH]) is satisfied, M 0 > N 0 , A is SOS if and only if vq > M 0 , and A is PSD 
if and only if Vq > N 0 . 

By Theorem Q1 we now only need to consider the case that rj(v 5 ,v e ) = 1. Suppose that 
for all v 5 and vq satisfying ((9]) , we have M 0 (v 2 , vq, v\, v 3 , v 3 ) = N 0 (v 2 , vq, vi,v 3 , V 5 ). Since the 
sets for PSD Hankel tensors and SOS Hankcl tensors are closed ra, this implies that for all 
v 5 and v 6 satisfying r](v 5 , v 6 ) = 1 , we also have M 0 (v 2 , v 6 , v 1 } v 3 , v 5 ) = N 0 (v 2 , v 6 , v 1 } v 3 , v 5 ), as 
long as iVo is defined there. Thus, in this case, by Theorem [3j there are no fourth order four 
dimensional PNS Hankel tensors under assumption (J2]). □ 

For the variables of M 0 and N 0 , we put v 2 and vq before v \, ^3 and v 3 , as v 2 ,vq play a 
more important role in the PSD and SOS properties of A, comparing with Vi,v 3 and n 5 . 

We now regard P = (v 2 ,v 6 , vi,v 3 , n 5 ) T as a point in !ft 5 . If M 0 (P ) = N 0 (P) } P is called 
a PNS-free point. We call 


S = {(v 2 , v 6 , vi, v 3 , v 5 ) t G 3? 5 : rj(v 5 , v 6 ) < 1} 

the effective domain. Theorem [5] says that if all the points in the effective domain are PNS- 
free, then there are no fourth order four dimensional PNS Hankel tensors with symmetric 
generating vectors. In the next sections, we will study more on PNS-free points. 

4 Theoretical Proofs of Some PNS-Free Regions 


4.1 Critical SOS Decomposition 

For the convenience, we present formally three ingredients used in theoretical proofs of 
this section. If a point belongs to the effective domain and enjoys these ingredients, it is 
PNS-free. 

Definition 1 Suppose that assumptions (0|) and (0|) hold and P = {y 2l vq , v \, v 3 , ns) T G 
S. Suppose that there is a number M such that A is SOS if v 0 = M, and a point x = 
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(xi, X 2 , x 3 , xfr 1 G 9ft 4 such that x\ + x\ > 0 and /o(x) = 0, where /o(x) = /(x) with vq = M. 
Then we call M the critical value of A at P, the SOS decomposition /o(x) the critical 
SOS decomposition of A at P, and x the critical minimizer of A at P. 

Theorem 6 Let P G S. Then P is PNS-free if A has a critical value M, a critical SOS 
decomposition /o(x) and a critical minimizer x at P. 

Proof Suppose that A has a critical value M, a critical SOS decomposition /o(x) and a 
critical minimizer x at P. Then we have M > MfrP) by the definition of M 0 . If vq < M , 
then 

/(x) = (u 0 - M){x\ + X4) + /o(x) < 0. 

This implies that NfrP) > M by the definition of N 0 . But N 0 (P) < M 0 (P). Thus, 
M 0 (P ) = N 0 (P) = M, i.e., P is PNS-free. □ 

We believe that all the effective domain S is PNS-free. In the next four subsections, we 
theoretically prove that some regions of S are PNS-free. 

4.2 A PNS-Free Segment 


Professor Man Kam Kwong pointed out that iV 0 (l, 1, 0, 0, 0) = 1, iVo(2,1,0,0, 0) = 8 
and i\T 0 (4, 0,0, 0,0) = 441, are integers. See also Table 1 in Section 6. He suggested us to 
considered these three points more carefully. Stimulated by Prof. Kwong’s comments, we 
derive the results of Subsections 4.2 and 4.3. 

We have the following theorem. 


Theorem 7 Suppose that P = (v 2 , u 6 , Vi, u 3 , u 5 ) T = (1,1, t, t, t) T , where t e [—1,1]. Then, 
P is PNS-free, with the critical value 1 and the critical minimizer (1,0, —1,0) T . 


Proof For P = (u 2 , u 6 , Vi, v :i , u 5 ) T = (1,1, t, t, t) T , where t G [—1,1], and M = 1, we have 


is SOS, and 


/o(x) = + x 2 + x 3 + x 4 ) 4 + —2~(x 1 ~x 2 + x 3 - x 4 ) 4 


/o(l, 0, —1, 0) = 0. 


Hence, P is PNS-free. 

By numerical experiments, we have the following conjecture. 


□ 


Conjecture 1 The segment L = {(v 2 , v 6 , v 3 , v 3 , v 5 ) r = (1,1 ,t,t,t) T : t G [—1,1]}, is the 
minimizer set of both M 0 and N 0 . 
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4.3 A PNS-free Planar Cone 


Theorem 8 Suppose that P = (v 2 , v e , Vi, v 3 , v 5 ) T = (v 2 , v 6 , 0, 0, 0) T with v 2 > v e > 1. Then, 

P is PNS-free. 

If we parameterize v 6 = b and v 2 = (9 +3b— 1)(0 2 + (3b — 2)6 — 3b+ 4). Then, the critical 
value at P is 

M =(9 +3b- T) 2 {39 2 + (106 - 6)9 + 36 2 - 10 b + 9) 
and the critical minimizer is x = (1, 0, —(9 + 36—1), 0) T . 

Proof Note that for v 2 > v 6 > 1, we may let v 6 — 6 and v 2 = (9+3b— l)(d 2 + (36—2)d—36+4), 
where the parameter 

9 > 9 = (6 — 1)5(6 + 1)1 + (6 — 1)§(6 + l)s — 26 + 1. 

In fact, 9 is the largest real root of the cubic equation v 2 — v 6 = 0. 

With the critical value as M = (9 + 3b— 1) 2 (3 9 2 + (106 — 6)9 + 36 2 — 106 + 9), the critical 
SOS decomposition at P is as follows 

/o(x) = —(v 0 xl + 2 v 2 XiX 3 + aixl ) 2 + ^-(v 0 xl + 2 v 2 x 2 X 4 + aixl ) 2 

v 0 n, 

+ ch 2 ((9 + 36 — \)x 3 x 3 + 3+)^ + tt 2 ((^ T 36 — 1 ) 3 + 3 + + X 2 ) 

6 1 \2 6(6 2 - 1) 2 

+ 7(3+3+ + X 3 X4 + 63+3+ + 63:1X4) H- - -(X1X2 + X3X4) 

6 6 

+ 6(u 2 - b)[x\x 2 2 + X3X4], 

where the involved parameters are as follows: 

a-! = —(9 2 + (46 — 2)9 + 36 2 — 46 + 1), 

2(9 2 + (46 — 2)9 + 6 2 — 46 + 4) 

° 2 “ 3 9 2 + (106 - 6)9 + 36 2 - 106 + 9' 

Since / 0 (1,0, — (0 + 36— 1), 0) = 0, the corresponding critical minimizer is x = (1, 0, — (# + 
36— 1),0) T . Hence, P = (i+, v 6 , 0, 0, 0) T with v 2 > v 6 > 1 is PNS-free. □ 

The cone C = {(v 2 ,vq,vi,v 3 ,v 3 ) t = (a, 6, 0,0, 0) T : a > 6 > 1} is a 45-degree planar 
closed convex cone. Its end point is just the mid point of the segment L = {(z+, v 6 , ++ 3 , ^s) 7 = 
(1,1 ,t,t,t) T : t G [—1,1]}, discussed in the last subsection. 

4.4 A PNS-Free Ray 
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In this subsection, we show that the ray R = {(u 2 , U6, ui, u 3 , U 5 ) T = (a,0,0,0,0) T : a < 
0} is PNS-free. Let a = —p, where p > 0 is a constant. We report that, at a point 
P = (—p, 0, 0, 0, 0) T , A has the critical value 

M = 3 \/ e 1 + 32v^+ — . 03 + 6 p 2 + 138p + 609, 

V 3^ + 32^ 

where 

9i := -p 6 + 272p 5 + 12608p 4 + 204032p 3 + 1558528p 2 + 5750784p + 8290304, 

0 2 : = —(p + 6) 2 (p + 4) 3 (p 2 + 4p — 16) 3 , 

0 3 ■= 9(p + 8)(p 3 + 152p 2 + 1728p + 5120). 

The function /o(x) enjoys a critical SOS decomposition: 

/o(x) = ^9fc(x), 

k =1 

where 

gi(x) = x 2 + 6x 2 x 4 + aix 2 + a 2 x 2 , 

g 2 (x) = x 2 + 6.T1X3 + a 2 x 2 + ccix 2 , 

g 3 (x) = a 3 x 2 x 4 + a 4 x 2 + a 5 x 2 , 
g 4 (x) = a 3 xix 3 + a 5 x 2 + a 4 x 2 , 
g 5 (x) = a 6 x 2 - a 6 x 2 . 


The involved parameters 

CX\ — 

are listed as follows: 

(p + 23)Mi(-p) - 9p 3 - 21p 2 + 105p + 9 

M 1 (—p) + 3p 2 + 6 p — 33 

«2 = 

-3 p, 

a 3 = 

\/—30 — 2 ai 5 , 

6(1 - 015 ) 

«4 = 

«33 

16p 

a 5 = 

«33 ’ 

«6 = 

/ 192p(a 1 5 — 1) 

V «33 


Theorem 9 Suppose that assumptions (HP and (HP hold. Then, for any constant p > 0, 
P = (—p, 0, 0, 0, 0) T is PNS-free. 
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Proof We only need to prove that there is a critical minimizer. Let 

X = (033, 0:35 + «36j — «35 — OL 36, — « 33 ) T - 

Then, we get <73 (x) = (x) = qs(x) = 0 immediately. Moreover, we have 

<?l(x) = q 2 (x) = (a35 + 0 36 ) 2 - 6(0:35 + « 36)«33 + « 15«33 “ 3 P a 33 = °- 

We check the validation of the last equality by a mathematical software Maple. Hence, 
/o(x) = 0 and x is a critical minimizer at P. Hence, we get the conclusion by Theorem [6] 

□ 

4.5 A PNS-Free Point 

We now show that the point A = (1, 0, 0, 0,0) T is PNS-free. In fact, the critical value at 
A is 

M = 477 + 3v/3906351 + 9120\/57 H -- - 74403 = . 

V 3906351 + 9120a/57 

The critical SOS decomposition of / 0 (x) is as follows 

7 

/o(x) = 5Z?fc(x) 2 , 

k =1 

= x 2 + 6x2X4 — 21x 2 + 04X4, 

= X2 + 6x4X3 — 21 x 2 + aqx 2 , 

= 2V3X2X4 + o 2 x 2 + 03X4, 

= 2 v / 3 xix 3 + o 2 x 2 + 03X4, 

04X4 0:4X4, 

= di X 1 X 2 + /3 2 ^1^4, 

= (3 4X3X4 + P2X1X4. 


where 

9i(x) 
92 (x) 
73 (x) 

94 (x) 

95 (x) 

96 (x) 

97 (x) 
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Some involved parameters are listed as follows: 


Pi 


P 2 


OL\ 


0-2 


Oi 3 


v /-6(M 2 - 36) (3M 2 - 4336) 
- 1302 M 2 + 25056 
fti(3ffi + 116) 

Pi + 12 ’ 

3 - 

22v / 3-^AA>, 

6 

8\/3 v+ 2 

Q “T o Pi’ 


0:4 — —42ai + 2 ck 2 g:3 + /?|. 

Theorem 10 Suppose that assumptions (Hj) and (E|) hold. Then, A = (1,0, 0, 0, 0) T is PNS- 
free. 


Proof Using the mathematical software Maple, we calculate 


/( x ) - 5^dfc(x) 

k =1 


-Pi - 120 pf + (4u 0 - 4944)/3 1 2 + 48n 0 - 69376 

4(A 2 + 12) 


(x\ + X 4 ). 


Substituting the value of v 0 = M and Pi, we get / 0 (x) — X]I=i 7fc( x ) = 0- 

Let x = (Pi,P 2 , — P21 — Pi) T ■ Obviously, we obtain g 5 (x) = g 6 (x) = <77 (x) = 0. We find 
that <73 (x) and <74 (x) vanishes if we rewrite all the parameters using p\. Using the value of 
each parameter, we hnd that q\ (x) = g 2 (x) = 0. Since x\ — Pi ~ 1.73, x is the critical 
minimizer. Therefore, this theorem is valid according to Theorem O □ 


5 Numerical Experiments 


We have proved in Section 4 that some regions are PNS-free. What about the other 
cases? We try to answer this problem by a numerical approach. We use the YALMIP 
software with an SOS module [15] [16] to compute M 0 (v 2 , v 6 , vi,v 3 , v 5 ), which is the smallest 
value of r>o such that the fourth order four dimensional Hankel tensor A with the generating 
vector (v 0 ,Vi,v 2 ,v 3 ,l,v 5 ,v e ,v 5 ,l,v 3 ,v 2 ,Vi,v 0 ) T is SOS. Gloptipoly [7] and SeDuMi [23] are 
employed to compute N 0 (v 2 ,v e ,vi,v 3 ,v 5 ), which is the smallest value of u 0 such that the 
Hankel tensor A is PSD. 
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V 2 \ v 6 

-.2 

-.1 

0 

.5 

1 

1.5 

2 

4 

-4.0 

3.54e4 

8.74e3 

3.76e3 

4.78e2 

3.12e2 

3.92e2 

6.23e2 

6.37e3 

-2.0 

2.98e4 

6.77e3 

2.73e3 

2.75e2 

1.25e2 

1.70e2 

3.57e2 

6.11e3 

-1.0 

2.72e4 

5.85e3 

2.26e3 

1.91e2 

6.15el 

9.26el 

2.73e2 

6.06e3 

-0.5 

2.59e4 

5.42e3 

2.04e3 

1.53e2 

3.78el 

6.41el 

2.48e2 

6.06e3 

0.0 

2.46e4 

4.99e3 

1.82e3 

1.20e2 

1.96el 

4.50el 

2.39e2 

6.07e3 

0.5 

2.34e4 

4.57e3 

1.62e3 

8.90el 

7.058 

4.18el 

2.45e2 

6.09e3 

1.0 

2.21e4 

4.17e3 

1.42e3 

6.21el 

1.000 

4.93el 

2.56e2 

6.11e3 

1.5 

2.09e4 

3.78e3 

1.23e3 

3.90el 

4.191 

5.69el 

2.67e2 

6.14e3 

2.0 

1.98e4 

3.41e3 

1.06e3 

2.02el 

8.00e0 

6.46el 

2.78e2 

6.16e3 

3.0 

1.75e4 

2.70e3 

7.28e2 

7.16e0 

1.66el 

8.01el 

3.01e2 

6.21e3 

4.0 

1.53e4 

2.04e3 

4.41e2 

1.23el 

2.60el 

9.60el 

3.23e2 

6.25e3 


Table 1: The values of M 0 (v 2 , v e , 0, 0, 0) = N 0 (v 2 , v e , 0, 0, 0) on some grid points. 

5.1 M 0 (v 2 , v 6 , 0,0,0) and N 0 (v 2 , v 6 , 0, 0, 0) 

First, we focus on two elements v 2 and v 6 of generating vectors and set V\ = v 3 = v 5 = 
0. By Theorem [H owing to the effective domain, we have b > — We choose v 2 = 
—4, —2, —1, —0.5, 0, 0.5,1,1.5, 2, 3,4 and = —0.2, —0.1, 0, 0.5,1,1.5, 2,4 and compute M 0 
and Nq in these grid points respectively. By our experiments, we found that these two 
functions are equivalent on all of the grid points. Thus, no PNS tensors are detected here. 
The detailed value of M 0 and N 0 are reported in Table Q] 

A more intuitional profile of M 0 = N 0 is illustrated in Figure [21 It is easy to see that 
(v 2 , Vq) = (1,1) is the minimizer of both M 0 and N 0 when we set V\ = v 3 = v 3 = 0. 

5.2 Nonzero Odd Elements of the Generating Vectors 

We consider the case that the generating vector of a fourth order four dimensional Hankel 
tensor has nonzero odd elements. According to Theorem 0 we say that v 3 and v fs must 
satisfy 17 ( 115 , vq) < 1. So we study them first and set v\ = v 2 = v 3 = 0. We compute a plenty 
of grid points with different v 3 and v§. The function M 0 (0, v§, 0, 0, v 3 ) is still equivalent to 
the function N 0 (0, v e , 0, 0, v 5 ). That is to say, no PNS tensors are found. 

The contour of M 0 (0, Ve, 0, 0, V 5 ) = M 0 (0, vq, 0, 0, V 5 ) is shown in Figure [3l We could see 
that the nonlinear contour of Mq = Nq = 500 looks like a fire balloon. 

Finally, we consider all of the elements of symmetric generating vectors of fourth order 
four dimensional Hankel tensors. The contours of M 0 (v 2 , v 6 , V\, v 3 , v 5 ) and N 0 (v 2 , v 6 , V\, v 3 , v 5 ) 
for various combinations of v 2 , Vq, V\, v 3 and v 3 are reported in Figure [H In all of our tests, 
values of the function M 0 (v 2 ,v e ,vi,v 3 ,v 5 ) in grid points are always equivalent to the corre- 
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Figure 2: The contour profile of M 0 (v 2 , v 6 , 0, 0, 0) = N 0 (v 2 , v e , 0,0, 0). 

sponding values of the function N 0 (v 2 , v^, V\, v 3 , v$). So, no fourth order four dimensional 
PNS Hankel tensors with symmetric generating vectors are detected. 

From Figures [3] and HI we could say that the second element V\ of the generating vector of 
a Hankel tensor affect functions Mq{v 2i Vq, V\, v 3 , V5) and N 0 (v 2 ,V 6 ,Vi,v 3 ,Vs) slightly. When 
we fix V 4 — 1, the middle element Vq of the generating vector v plays a more important role 
since it has direct impact on the effective domain. 

6 Final Remarks 


In this paper, we investigated the problem whether there exist fourth order four dimen¬ 
sional PNS Hankel tensors with symmetric generating vectors. Theoretically, we proved 
that such PNS Hankel tensors do not exist on the segment L = {(v 2 ,v e ,Vi,v 3 l v 5 ) T = 
(1,1 ,£,£,£) T : t E [—1,1]}, the cone C = {(v 2 ,v^,Vi,v 3 ,v 5 ) t = (a,b, 0,0, 0) T : a > b > 1}, 
the ray R = {(v 2 ,vg,vi,v 3 : V 5 ) t = (a,0,0,0,0) T : a < 0} and the point A = (1, 0, 0,0, 0) T . 
The critical value on L is simply 1. It is interesting to note that the critical values on C are 
a polynomial in an auxiliary parameter 9 with degree four. However, the critical values on 
R and A are irrational. This indicate that a complete proof that fourth order four dimen¬ 
sional PNS Hankel tensors with symmetric generating vectors do not exist may not be easy. 
However, numerical tests also indicate that such PNS Hankel tensors do not exist. Thus, we 
believe that there are no fourth order four dimensional PNS Hankel tensors with symmetric 
generating vectors. 
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Figure 3: The contour profile of M o (0, v 6 , 0,0, v 5 ). 
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